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Abstract 



Critical examination is made on the relation between the mass shift of 
vector mesons in nuclear medium and the vector-meson — nucleon scattering 
length. We give detailed comparison between the QCD sum rule approach by 
two of the present authors (Phys. Rev. C46 (1992) R34) and the scattering- 
length approach by Koike (Phys. Rev. C51 (1995) 1488). It is shown that 
the latter approach is mortally flawed both technically and conceptually. 
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I. INTRODUCTION 



The question of changes in the vector meson properties in nuclear medium is of interest 
in relation to the physics of nuclear matter |l] and relativist ic heavy ion collision 0. In 
particular, if there are spectral changes of vector mesons in medium, it could be observed 
directly through the lepton-pair spectrum in j—A, p—A and A— A reactions ||. The problem 
has been studied in effective hadronic models and QCD sum rules (QSR) generalized to finite 
baryon density. In QSR, it was shown that the vector meson mass drops to about 10—20 
% of its vacuum value in nuclear matter density Q. This is consistent with the idea of 
the Brown-Rho scaling || and also with Walecka model calculations including the vacuum 
polarization ||. 

Recently, Koike J?| claimed that the previous QSR calculation for vector meson in nuclear 
medium j|| is not correct and that the vector meson mass should increase slightly. Since 
both ref. 0] and M on the vector meson masses are based on the Fermi gas approximation 
for nuclear matter, it is desirable to clarify the origin of the difference. In this work, we 
will first give a correct account of the Fermi-gas approximation of the QSR in medium and 
show dangers to relate mass shifts with the scattering length. Secondly, we will show that 
the approach in M is technically wrong so that in a correct application of QSR in medium 
Q the vector meson mass does decrease in nuclear medium. 

This paper is organized as follows. In section II, we will summarize the results of ref. Q 
and make new remarks which are not given in jij but are relevant for the later discussions. In 
section III, the essential points of the Koike's claim are summarized. Section VI is devoted 
to critical investigation of p]]. We will show that his method is mortally flawed. In section 
V, further elaboration of ref. || is discussed. 



II. QSR IN NUCLEAR MEDIUM 

Let's start with the retarded current correlation function in nuclear matter, 

II* (w, q; n B )=ij dV^RJ^x) J v (0)) n , (1) 

where g M = (u>, q) and RJ At (x) J u (0) = 9(x°)[J^(x), J u (0)] with the source currents defined 
as = \(u'y fl u^fd'y fl d) (— (+) is for the p (co>)-meson). n B denotes the baryon number density 
and (■)„ is the expectation value in the ground state of nuclear matter. 

Although there are two independent invariants in medium (transverse and longitudinal 
polarization), they coincide in the limit q — ► and reduce to II* / (— 3uj 2 ) = IT*. IT* satisfies 
the following dispersion relation, 

^ T-r P / 9x 1^ f°° , nlmH R (u) , , . . 

ReU R (u 2 ) = -P / du 2 — V + (subtraction). (2) 

n Jo u 2 — ijj 1 

In QSR, the spectral density Imll^ is modeled with several phenomenological parameters, 
while Rell^ is calculated using the operator product expansion (OPE). The phenomenolog- 
ical parameters are then extracted by matching the left and right hand side of (fj) in the 
asymptotic region uj 2 — > — oo. 
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Although the nuclear matter ground state has complicated structure, the Fermi-gas ap- 
proximation taking into account the Pauli principle among the nucleons is a good starting 
point. In this approximation, Er^ reads 

/Pf d 3 v 
— — ^-T^(a;,qip), (3) 

with 

T*(u,q | p)=ij rf 4 xe^(iV(p)|RJ At (x)J,(0)|iV(p)). (4) 

Here 7 is a degeneracy factor (7=4 in nuclear matter), |iV(p)) is the nucleon state with four 
momentum p^ = (En,p) in the covariant normalization (iV(p)|JV(p')) = (2tt) 3 2En5 3 (p — 
p'). The spin average for the nucleon state is assumed in (Q). T^(u,q | p) is nothing but 
the spin-averaged forward scattering amplitude of the external current J M and the nucleon. 
In nuclear matter, p is integrated out in the range 0<|p|<p/(pyis the fermi momentum). 

The OPE for ReU R (uj 2 ), which is the same with the OPE for the causal (Feynman) 
correlation U f (uj 2 ), has a general form at u 2 = — Q 2 — > —00, 

Ren > 2 ^-00)=E^W 2 ^ 2 )(^(^))n , (5) 

where \i is the renormalization point of the local operators, which separates the hard scale 
\uj\ and soft scales such as Aqcd and pf. The Wilson coefficients do not depend on the 
medium effect and only the expectation values {Oi(fj, 2 )) n have the pf dependence. 

The local operators Oi(fi 2 ) in the vector meson sum rule are essentially the same with 
those in the lepton-nucleon deep inelastic scattering (DIS) and can be characterized by their 
canonical dimension (d) and the twist (r=dimension-spin). They are given in |IJ up to 



dimension 6 operators and we will not recapitulate them here. Since we are taking q — > 0, 
eq. (|J) is an asymptotic series in 1/uj 2 or equivalently an expansion with respect to d. 

Now let us examine the density dependence of (Oj(/i 2 ))„. In the Fermi-gas approxima- 
tion, it reads 

(0 >»=< 0> ° + ^ P '(sfk (Ar(p)|0|jV(p)) ' (6) 

For the scalar matrix elements such as (N\qq\N) and (N\G 2 \N), one can utilize the 
information of the tt — N sigma term E^jv an d the QCD trace anomaly |TTJ . The results are 



Mn = (uu) + ■ n B ■ B x (x), (—G 2 ) n =(—G 2 ) --mP-n B -B 1 (x). (7) 
2m 71 7r 9 

Here m(l GeV)= (7 ± 2) MeV is the average value of the current quark masses of u and d 
quarks. The parameters we use are T, nN = (45 ± 7) MeV, (^-G 2 ) = (350MeV) 4 , (uu) = 

(-230MeV) 3 § and mf = 770 - 830 MeV with y = 2(N\ss\ N)/((N\uu\ N) + (N\dd\N)) = 
0.22 — 0.12 PJTT[]. Non-leading pj corrections are contained in Bi(x) defined as 

B x {x) = ^-(xVl + x 2 - ln[x + Vl + x 2 ]) ~ 1 - ^-x 2 + ^-x* ■ ••, (8) 
2x 6 10 00 
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where x = pf/m^ and n B = 7p 3 /(67r 2 ). For nuclear matter density uq = 0.17/ fm 3 (x = 
0.27), x dependence of Bi(x) can be safely neglected. Note also that (dd) n = {uu) n , since 
we are considering the symmetric nuclear matter (N=Z). 

As for the four quark condensate in medium, there are no experimental data available 
yet. Thus we use a simple mean field approximation in nuclear matter Jlj; 

((qi^ a q) 2 )n * -m,\ a q) 2 )n * y [{(qqf)o + 2(qq)o(N\qq\N) ■ n B • B 1 (x)} , (9) 

As for ({qq) 2 )o at 1 GeV scale, we will use the canonical value (— 281MeV) 6 [|L5] with 
a s (lGeV) ~ 0.36. This number, which is substantially larger than the current algebra 
value (qq)o — (—230 MeV) 3 at 1 GeV, should be considered as an effective one containing 
non- leading 1/N C contributions. As for (qq)o in the second term of @, it is not clear whether 
one should use (— 230MeV) 3 or (— 281MeV) 3 . (Note that the latter number was used in Q.) 
Taking into account such ambiguity as well as the "experimental" errors of E^jy and rh, we 
adopt (qq) (N \ qq \ N) = (-256MeV) 3 • (45/14) • (1 ± 0.368) as a standard value to be used 
in (|). 

The twist- 2 quark bilinear operators with dimension 4 and 6 have non- vanishing matrix 
elements in the medium. Their nucleon matrix elements are related to the parton distribution 
function in DIS as, 

(N^STq^D^ ■ ■ ■ D, n q\N(p))l 2 = (-i)- 1 ^T,,^ 

AlJn 2 ) = 2 f 1 dxx n - l [q(x, /j, 2 ) + q(x, fi 2 )]. (10) 
Jo 

Here T Atl ... /Jn = [p w • • -p^ n — (trace terms)]/2mAr and ST makes the operators symmetric and 
traceless. For the parton distribution function q(x, /i 2 ), we take the LO scheme in ref. [12] at 



the scale fi 2 = 1 GeV which is close to the relevant Borel mass. Then we obtain A\ +d ~ 0.9, 
j^u+d ^ Q.12 at 1 GeV. The fermi-motion corrections are contained in B^ix) = \/l + x 2 and 
B 3 (x) = yi + x 2 {l + |x 2 ) for spin 2 and spin 4 operators respectively. The deviation of 
B 2 ^{x) from 1 is again small at nuclear matter density. 

In reference [§, all the operators up to d=6 except for relatively small twist-4 spin-2 
operators are taken into account. In section IV, we will discuss the effect of the twist-4 
operators to the result of 0] . 

In the vacuum QSR, the spectral function (i.e. Imll^ in eq.(fj)) is modeled with a 
resonance pole and the continuum. In the medium, we have to add additional singularities 
below the lowest resonance pole within the Fermi-gas approximation, which is called the 
Landau damping contribution ||13|| . For q — > 0, it is calculable exactly and behaves like a 
pole at uj 2 = (see Appendix A for the proof). In total, the hadronic spectral function looks 
as 

8-KlmIl R (u > 0~) = 5(u 2 ) Psc + F5(u 2 - m\) + (1 + — )9{u 2 - S ) (11) 

7T 

= Phad.iu 2 ), 



with p sc = 2n n B / ^p 2 + m 2 N ~ 2n ns/mN- rny, F and So are the three phenomenological 
parameters to be determined by the sum rules. 
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Matching the OPE side and the phenomenological side via the dispersion relation in 
the asymptotic region uj 2 — > — oo, we can relate the resonance parameters to the density 
dependent condensates. There are two major procedure for this matching, namely the finite 
energy sum rules (FESR) |14| and the Borel sum rules (BSR) JIBf which are summarized as 
the following forms: 

/ dsW(s) [p had .(s)-p OPE (s)]=0, (12) 
Jo 



W{s) 



s n 6(S -s) (FESR), 
(BSR). 



s/M 2 



Here the spectral function Phad.( s ) stands for eq. ([□]). p OPE (s) is a hypothetical imaginary 
part of 11^ which, through the dispersion relation (|2|), reproduces eq.(^). For more details 
on (IP) and the explicit form of p OPE , see section 2 of ref. fT0|j. 



FESR and BSR for U R {u 2 ) 

First, for the qualitative argument, let us write down the FESR for the rho (omega) 
meson in the chiral limit. This can be easily obtained by taking the first three moments 
n = 0,1,2 in (Pj): 



S (l + —) = -2n 2 m- 1 -n B (n = 0), 

7T 



Fm 2 - (1 + —) = -Q 4 - 2K 2 A u 1 +d m N n B = -Q 4 (n = 1] 



o2 

( 




= -Qa 


1 + -) 


2 1 


7T 




^0 / 




= -Qe 




3 1 


7T 




(fG 2 ) n 





4 - ^(1 + ^) = -Qa " yvr 2 A^> B = -Qe (n = 2), (13) 

where Q 4 (Q 6 ) is 2 ^(^G 2 ) n (^7r 3 (a s (gg) 2 ) n ) taken up to linear in n B . Using the three 
relations above, we can determine the three phenomenological parameters F, So and my 
or equivalently, the changes from the vacuum values SF, 5 So and 8m v . Important density 
dependence comes from A^ +d n B and the the 4-quark condensate Qq, which can be shown 
by solving eq. ([13D numerically. 

Although both FESR and BSR give the same qualitative result, BSR is more reliable for 
the quantitative estimate since it is rather insensitive to the assumption on the continuum. 
The rho (omega) meson mass in the BSR is given as 



m 2 



_ (1 + g ) (1 - e- So/M \l + m) - + ikQa 

M2 (l + t)(l-e-^ 2 ) + ^Q 4 -^Q 6 -p sc ' 1 J 

In Fig.l, the Borel curve (m v — M 2 relation) is shown for different baryon densities. The 
continuum threshold So is chosen to make the Borel curve as flat as possible in the Borel 
window M^ in < M 2 < M^ ax at given density. We take density-independent window M^ in = 
0.41GeV 2 and M 2 max = 1.30GeV 2 in our analyses. More general procedure with density- 
dependent window (see e.g. section 4 of ref. [TO]) does not change the results quantitatively. 
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Fig.l 

By making a linear fit using the values at n = n Q and n = n B , we get 



m v (0) 



So(n B ) 
\ S (0) 
F{n B ) 
F(0) 



(0.16 ±0.06) — , 

n 



(0.15 ±0.05)^ 
n 



(0.24 ±0.07)-^ 

n 



(15) 
(16) 
(17) 



These values are slightly different from our previous ones in 0] where the uncertainty dis- 
cussed below @ is not taken into account. For the decreasing rho (omega) mass, the twist 
2 and the scalar matrix elements are equally important. 



III. MASS SHIFT AND THE SCATTERING LENGTH 

In this section, we will first summarize the claims given explicitly or implicitly in J7|. 
Although they look plausible in the first look, every statement is invalid as we will show in 
section III.A-C. 

(A) Koike starts with eq.@ and makes a low density approximation for the second term 
fPf d 3 p , T*(u;,q|0) , , 

This corresponds to the assumption that all the nucleons in nuclear medium are at rest 
(p = 0). If one further takes the kinematics q = and uj ~ m v , T R = (u, q = | 0) is 
written by the V — N scattering length as 

r «^3F^ 24,(m +m K, + 
8^^ [uj z — m l v Y 

where a VN = (a x / 2 + 2a 3 / 2 )/3 with au 2 and a 3 / 2 being the V — N scattering length in the 
spin-1/2 and spin-3/2 channel respectively. 1Z{uj 2 ) is the term less singular than the leading 
double-pole term in the Laurent expansion around u 2 = rriy. 

By substituting these expressions to eq.(0) and take the leading term in eq.(^), one 
arrives at the formula 

1 t7 r i n \ Fm 2 v ( 1 l2na VN (m N +m v )/m N \ 
^>^m l ,,q = 0;?i ij ~ — -f — - + — — n B , 

oc — 1 — (20) 

uo z — [m v + om v Y 

with 
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m N +m v 

dm v = 6vr a VN ■ n B . (21) 



m N m v 



Hence the positive (negative) scattering length gives an increasing (decreasing) mass in the 
medium. 

(B) To estimate the magnitude and sign of the scattering length a VN in fl2"ip, Koike formu- 
lated QCD sum rules for T R . He starts with the unsubtracted dispersion relation; 

Re T r {uj 2 ) = - / du 2 — ^. (22) 



IT JO V? — UJ 2 



The OPE for T R is expanded up to 0{1/Q A ) as 



where cx^/2m N = dQ^ 6 /dn B at n s =0. The absence of the logarithmic term in ( p3|) indicates 
that subtraction is not necessary in fl22j) . Motivated by eq . (|T9|) , the imaginary part is 



parameterized as 

8nlmT R (u > 0) = M'(« 2 - m v (0) 2 ) + 6 2 5(u 2 - m v .(0) 2 ) + b 3 5(u 2 - S Q (0)), (24) 

with three unknowns 61,2,3 an d known vacuum parameters my(0) and So(0). (Note that 
61 = Fm 2 v 2An{m N + m v )a VN .) Above parametrization is equivalent to taking the following 
ansatz for the real part 

BeT R (O 2 ) - bl b 2 b 3 

Ke ±hadW } (m v (0) 2 + Q 2 ) 2 m v (0) 2 + Q 2 + S (0) + Q 2 ' 1 } 



By constructing a borel sum rule using (0)— (|24|) , Koike obtains positive scattering length 
a p (aj) ~ 0.14 (0.11) fm, from which he concluded that the mass shift 5m v in eq.([H|) must be 
positive. 

(C) Koike further claims that the above procedure (A) + (B) is equivalent with doing the 
medium sum rule for but not for n^. (Note that 11^ = —3u) 2 U R when q = 0). This 
can be seen as follows. The dispersion relation in medium for 11^ reads 

ReU R (cu 2 ) = -P I" du 2 lm[U ^ ] + (subtraction). (26) 
w ir io u 2 — uj 2 

If one adopts eq. ([TT|) for Imn^ and uses the relation H R = —3uj 2 U r , one obtains 

1 a 
(--)87rIm[H* («)] = Fm 2 v 5(u 2 - m 2 v ) + (1 + ^)u 2 6{u 2 - S ), 

6 7T 

= u 2 p had .(u 2 ). (27) 
Since u 2 5{u 2 ) = 0, the Landau damping term in ( |TT1) does not arise in (f2~7l). 
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Expansion of the l.h.s. of ( p6|) in terms of % gives [ub independent term] + [n^x 



eq. (p3|) 1, while the same expansion of (|27| ) gives [% independent term] + [n B x eq.(f24|)]. The 
latter is obtained simply by writing m v = m v (0) + 5m v , S = S (0) + 5S , F = F(0) + SF, 
expanding ( f27| ) up to linear in 5m v , 6F and 5 So and doing the following identification. 



b x j2m N = -Fm 2 v 5m 2 v , b 2 /2m N = m 2 v 5F + F5m 2 v , b 3 /2m N = -S 5S . (28) 

This means that the sum rule for T R eq.(|2"2"|) is equivalent to the linear density part of 
the sum rule for 11^ eq . (|26"D . Assuming that his procedure (A) + (B) is right, Koike thus 
concluded that (i) the medium QSR using 11^ must give the increasing vector-meson mass, 
and (ii) the result of the medium sum rule using U R in || must be wrong. 

In the following, we will discuss that each of the above arguments is invalid. The sub- 
section numbers III. A, III.B and III.C in the following correspond to the statements (A), 
(B) and (C) in the above, respectively. 



A. Mass shift and the scattering length 

First of all, let us show that eq.flTED is a useless approximation around nuclear matter 



density. In eq.([L8|), the motion of nucleons and the Pauli exclusion principle in nuclear 
matter are completely neglected. Such approximation is valid only when (i) the nucleon 
density is extremely low, or (ii) T R u (uo,(\ = | p) is almost constant as a function of p in 
the interval <| p | < pf. 

Since we are not interested in the case (i), let us concentrate on (ii) and see whether (ii) 
is plausible or not. At nuclear matter density, the fermi momentum is sizable Pf — 270 MeV. 
Thus we should consider e.g. the p — N scattering from a/s = m p + m^ = 1709 MeV through 

y/s = {{m p + ^m 2 N + p 2 ) 2 — p 2 ) 1 / 2 = 1726 MeV. In this interval, there are at least two s- 
channel resonances A(1710) and A(1720) and also there are two nearby resonances just 
below the threshold A(1700) and A(1700) fIB|| . They are all possible to have coupling with 
pN. This means that T R V has a rapid variation as a function of | p | between | p | = and 
| p |= pf due to the effect of these s-channel resonances and it is impossible to approximate 
it by the threshold value (i.e. the V — N scattering length) T r (uj = my, q = | p = 0). See 



Fig. 2 for a schematic illustration of the s— channel contributions [16 



Fig. 2 



What one can expect at best is the approximate linear density formula written in terms 
of the average of T R U in the region <| p \< Pf- 

[Pfd*p R / T R (u, q =0\p) \ 

y (27j^ T - ( "' q=0|p) " nB y^e- n — /' (29) 

where (■) stands for the average over the above momentum interval. 

It is easy to see what is wrong in ( |I8D in physical terms: the rho-meson at rest (q = 0) in 
nuclear matter will suffer the scattering from the nucleons having various three-momentum 
p in the interval <| p |< pf. The net effect should be the averaged strength of the 
scattering and not the scattering length defined at p = 0. This is particularly so when T R U 
has a rapid p dependence. 



S 



The above point is well known for the nucleon in nuclear medium. The optical potential 
for the nucleon at rest in nuclear matter cannot be approximated by the N-N scattering 
length multiplied by the nuclear density. In fact, the N-N forward scattering amplitude 
T NN (p) has a huge momentum dependence due to the the deuteron state and the "almost" 
bound state near the threshold. The relevant quantity for the nucleon optical potential is 
not the scattering length but the averaged scattering amplitude in the interval <| p \< Pf- 
This point was recently emphasized by Furnstahl and one of the present authors M: They 



criticise the paper by Kondo and Morimatsu [IT] who use a similar approximation with (|18|) 
to analyse the nucleon in nuclear medium. 

It is now clear that the mass shift and the scattering length cannot be related directly at 
nuclear matter density. In the approach of ref. |§], neither eq.([0|) nor eq. (|29"D are adopted, 
thus one does not suffer from this problem. 



B. Is scattering length calculable in QSR? 

Since eq.(|l"8"l) and hence eq . ([2lD are not valid around nuclear matter density, it is useless 
to relate the scattering length with the mass shift. Nevertheless, the V — N scattering 
length itself could be an interesting physical quantity to be calculated in QSR. We will 
show, however, that the method in ref. J7| to estimate a VN is erroneous. 

Let us look at (|23|) and ([25|) and compare them at Q 2 — > oo, which corresponds to 
the FESR for T R . One immediately realizes that only two independent equations can be 
obtained 

b 2 + h = ci, 

6 1 -m l/ (0) 2 6 2 -S , (0)6 3 = -c 2 , (30) 

whereas one needs three equations to solve 61,2,3- This happens because the OPE is calculated 
only up to 0(1/ Q 4 ) in (^3|). ( p0|) clearly shows that it is impossible to predict the scattering 
length a VN (which is proportional to 61). 

In QCD sum rules, one should always check that the number of phenomenological pa- 
rameters to be determined is equal or smaller than the number of OPE terms, otherwise 
sum rules are not closed and cannot give predictions. The FESR provides a useful tool to 
do this consistency check. 

If one tries to make Borel analyses without the consistency as has been done in ref. [|7| , 
one simply obtains a fake result and does not get any stability of the Borel curve. By blindly 
doing a Borel sum rule without looking for a stability region, Koike obtains a positive number 
for b\, in which it is implicitly assumed that the contribution from dimension 8 operators is 
zero. In terms of finite energy sum rule, this amounts to the following condition, 

2m v (0) 2 b x - m v ,(0) 4 b 2 + S (0) 2 b 3 = 0. (31) 

With this assumption and eq. (|3~0"D , we have three equations and the unknown constants can 
be determined. Using the fact that c\ ~ 0, c 2 > and So > m 2 v) one finds a positive value 
for b\. However, the assumption eq.(^TJ) has no ground and we do expect non-negligible 
contribution from the dimension 8 operators. Since it is technically very hard to calculate 
these dim. 8 condensates ( l/Q 6 terms of T R ), it is almost hopeless to get reliable a VN in 
QSR. 
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C. Use of QSR for II* 



Here, we will explicitly demonstrate that the sum rules for n* recommended in |7j does 
not work at all. Let's first start with the FESR for n* ■ 

jf* ds s n pmll* (a)^ - ImnJ M ( S ) OP£ ] =0 (n = 0, 1, 2, • • •)• (32) 

Then one immediately finds that only two relations corresponding to n — 0, 1 are obtained, 
and they turn out to be equivalent with the second and third relations in eq. (|T3"D . (Note that 
n=2 in (^) cannot gives another condition because OPE is calculated only up to dimension 
6 operators.) 

(l + -) = -Q 4 , (n = 0) 

(1 + -) = -& (n = l). (33) 

3 7T 

There are three unknowns F, m v , So, while only two relations are available. Thus unless one 
introduces extra assumption, it is impossible to solve for three even in the vacuum. This is 
exactly the same problem which we have discussed in sec. III. B. 

One should also note that the missing condition is the duality relation for the spectral 
density: 

So 

ds [phad.(s) - p OPE (s)} = 0, (34) 

with phad. and p OPE being the spectral densities for n*. This local duality is the cornerstone 
of vacuum QSR and holds also in medium since there are no dimension 2 operators in OPE. 
If one wants to get reliable result from FESR, One has either to work out the OPE up to 
dimension 8 operators (which is a formidable task) or to start with n* as in ref. M. 

Here one may ask that "why not take n = — 1 moment in (^) to obtain another re- 
lation?". Such procedure, however, introduces an ambiguity at s — 0, since one can add 
any function proportional to s5(s) to Imn* which does not modify n — 0, 1 sum rules but 
modifies n — — 1 sum rule. If one tries to remove this ambiguity, it is necessary to start 
with Ii R and to evaluate the Landau damping term in (|ll|) which is exactly calculable as we 



o 



mentioned ||I9 



Let us now turn to the BSR for n^„: 

ml 2(1 + SO fl - e -WM 2 (i + + Jk,)) _ 1q 



V 



M2 (1 + t) 1 - e~ s °/ M2 (l + M - + 



(35) 



Eq.(|l4D and eq . (p5|) in this paper correspond to eq.(13) and eq.(14) in respectively. Since 



Q 6 enters with opposite sign in the r.h.s. of eq.(|T4D and eq.(|33D and Q 6 decreases in medium, 



Koike simply concluded that m v decreases in (0), while it increases in fl3"5|) . This conclusion 



is too naive: in fact, Sq in the r.h.s. of these equations is also density dependent, which can 
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change such naive expectation, density dependence of So can be in principle determined by 
the Borel stability procedure we have discussed in section II. 

In order to see whether this procedure works or not for eq. fl35|) , we have shown the Borel 
curves for different values of So in Fig- 3(a) (at zero density) and in Fig. 3 (b) (at nuclear 
matter density). Fig. 3(a) shows that the Borel curve does not have any plateau in the 
relevant range of M 2 (say 0.41 < M 2 < 1.30), which implies that one cannot determine So 
and hence m v even in the vacuum. The situation is the same at the nuclear matter density 
as is shown in Fig. 3(b). Again, one cannot determine So (no) by the Borel stability method 
and hence m v (no), which implies that there is no hope to determine the mass shift at finite 
density. If one sticks to a specific value of So (say 2.0GeV 2 ) and uses it at any density, one 
finds that m v (n B )/m v (0) > 1 for given M 2 from Fig.3. This is equivalent with the "naive" 
(and wrong) argument by Koike. The correct procedure is to compare m v {n B ) (calculated 
with So(n fl )) and m v (0) (calculated with So(0)). 

The "bad" Borel curves in Fig. 3 is quite in contrast to the "good" Borel curves for 11^ 
in Fig.l. The latter shows beautiful stability in the vacuum as well as in the medium, which 
makes one possible to determine So(n B ) at each density and hence m v (n B ). 

The reason of the failure of the BSR for H R is twofold. Firstly, the higher dimensional 
operators in OPE is rather important for sum rules. We have already seen this in (|33"D 
where n = 2 sum rule can be obtained only when one has dim. 8 operators in OPE. In BSR, 
the lack of the information of dim. 8 operators arises as a instability of the Borel curve at low 
M region. Inclusion of the dim. 8 operator would make the curve more flat. Secondly, the 
continuum contribution is more important in than in U R , since the spectral function is 
increasing linearly in the former case. This makes the prediction of the resonance parameters 
less reliable in the former. 

Let us summarize here the lessons we learned in subsections III. A, III.B and III.C. Firstly, 
the mass shift and the scattering length does not have direct relation in nuclear matter due 
to the momentum dependence of the V — N forward scattering amplitude. Secondly, sum 
rules for the V — N scattering amplitude cannot predict the V — N scattering length without 
dimension 8 operators in OPE. Thirdly, sum rules for U R does not work at all even in the 
vacuum without dimension 8 operators in OPE. Thus all the claims in ref. J7| are shown 
to be erroneous. Also, only the consistent QSR in medium currently available is the one 
starting from U R given in ref. ||. 



IV. SEVERAL COMMENTS 

Full twist 4 calculation 

There are three kinds of twist-4 spin-2 operators contributing to the p, u sum rules. At 
present, their nucleon matrix elements are not known. However, the values of two different 
combinations in the transverse and longitudinal structure functions of the nucleon have been 
obtained by two of us pOl^Tf by analyzing the recent DIS data at CERN and SLAC. Let us 



further make the following assumption 
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(dT^A u d) _ (dj^Dyd) _ 
(uT^&vu) {u^^Dyu) 

where V ^ is some gamma matrix and A M an isospin singlet operator. Then, it is possible to 
uniquely determine the nucleon matrix element of twist-4 spin-2 operators appearing in the 
p, u sum rule from combinations of experimental values. It gives the following contribution 
to eq.flD, 

QlP) = ~ [iq ' P ^ mW] (-(! + ml + {ft + \K) + K ud ± K ud ) (37) 



The K l u i = 1,2, g are defined in ref. |2D| and — (+) corresponds to the u>(p) case. 
Choosing (3 = 0.476 as in ref. [ 20 1 the value inside the large round bracket is A 4 = 0.40 



(0.24) GeV 2 for the u (p) meson. 

Now the effect of twist-4 matrix element can be estimated by making the following 
substitution for the dimension six operators. 



;Q 6 + ^n 2 A u 3 +d m 3 N n B ) -> (Qe + fn 2 A u 3 + d m 3 N n B + 2n 2 m N A*n B ) (38) 



The net effect of twist- 2 + twist-4 is estimated to be 2.36 (3.29) times larger than the 
twist 2 effect alone in the rho (omega) channel. This could change the slope of the mass 
shift in ( |I5D from 0.16 to 0.10 (0.075) for the rho (omega) meson. Further investigation is 
necessary however to draw definite conclusion on the magnitude of the twist 4 effect. 

Fermi momentum correction 

The small higher density effects coming from fermi-momentum correction can be es- 
timated by looking at b{(x). One easily finds that Si (0.27) = 0.979, S 2 (0.27) = 
1.036, S 3 (0.27) = 1.157 and Si (0.34) = 0.967, S 2 (0.34) = 1.056, S 3 (0.34) = 1.251, where 
x = 0.27(0.34) corresponds to nuclear matter (twice of nuclear matter) density. Thus the 
effect can be safely neglected at nuclear matter density. 

Possible new structure in Imll^ 

It is possible that the density dependent change of the OPE side is balanced by some 
new structure appearing in the spectral density below the resonance mass. In QCD sum 
rule approach, it has to be included by hand before matching to the OPE. Such possibility 
has been examined by Asakawa and Ko |22| by redoing the medium QCD sum rules for 



the vector meson including other complex structure of the spectral density in the nuclear 
medium induced by the tt, A, N, p dynamics. They found that even in that case the vector 
meson mass has to decrease in order to be consistent with the OPE side. 



4> meson sum rule 

The formalism for calculating the change of <fi meson sum rule is the same as that of 
the p and u>. However, in the (f) case, one must include the effect of the strange quark mass 
in the OPE and this will introduce some basic difference || . In the p, uj sum rule, the 
density dependence in OPE is dominated by (d, r) = (4, 2) and (d, r) = (6, 0) operators. 
However, in the case of the <fi meson, the dimension 4 strange quark condensate (m s ss) is 
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not suppressed either by 1/ 4ir 2 or the light current quark mass and consequently dominates 
the OPE. In the medium, the change of this condensate, which comes from the the K-N 
sigma term, dominates the small changes in other condensate and introduce a non-negligible 
mass reduction by 3 ~ 5% J| . It is amusing to compare this result with that of an effective 
model calculations |23| in which the K-N sigma term also induces a small reduction of the 
6 mass. 



V. SUMMARY 

To estimate the mass shift of vector mesons in medium, we have carried out a detailed 
comparision between the approach based on the modification of the vacuum QSR [|J and 
that based on the scattering length . We have shown that the latter approach is erroneous 
by the following reasons: 

(i) The mass shift and the scattering length does not have direct relation in nuclear matter 
due to the momentum dependence of the V — N forward scattering amplitude. 

(ii) Sum rules for the V — N scattering amplitude cannot predict the V — N scattering length 
if dimension 8 operators in OPE are not included. 

(iii) Sum rules for IPj have no predictive power both in the vacuum and in the medium if 
dimension 8 operators are not included. 
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APPENDIX: A 



Here, we derive the Landau damping term p sc . In the Fermi gas approximation, the 
spectral density has two types of contribution. The annihilation term, which is non-zero 
above the two particle threshold uj 2 > q 2 + 4m^, and the scattering term, which is non-zero 
in the space-like region uj 2 < q 2 . We are interested in the second term. For finite q, the 
spectral density, contributing to the longitudinal polarization, can be obtained by looking 
at the 00 component of the imaginary part of eq. ([]]). 



Imrioo q->o Imn w 

/^,q) = — ^ — ^-r- (Ai) 



Looking at the spectral representation, it is easy to identify the following scattering 
contribution, 

(2tt) 4 f d 3 k x d 3 k 2 2 



"•^ =~Wl (2,)3 gl (2,)3k l<JV(k ' )|J ° |JV(k2))l (A2) 



xS(uj -E x + E 2 )S 3 (q - kx + k 2 ) [n F (E 2 ) - n F (^i)] . 



Here, Ei = y k 2 + m 2 N (i=l,2), and the np(Ei) = 0(Jk 2 + m 2 N — Ei), where kf is the fermi 
momentum. 

In general, the nucleon expectation of the isospin current has two form factors. 

(iV(k 1 )|j;|iV(k 2 )) = fiflcOy [FMl, + F 2 (q)ia^\ u(k 2 ). (A3) 



Substituting this into eq . (|A2|) , one obtains 

7 



pi \u, qj 



256tt 2 



/ dx \{l-x 2 )E?{q)+x 2 {u 2 -cf)F 2 {q)} (A4) 
x ( 2 „ F ( WE+m) - 2n F ( ) Jtf - u,*), 



where v = [1 — 4m 2 / (oo 2 — q 2 )] 1 ^ 2 . 

Now, we want to take the limit |q| — > 0. In this limit, the constraint < uj 2 < q 2 
also forces uj to approach zero. Consequently, the contribution proportional to F 2 vanishes, 
because it is multiplied by either q 2 or uj 2 . As for the other term proportional to Fi(q), the 
integral becomes increasingly large as |q| — > such that the integrated quantity of pf(co, q) 
within the phase space for uj remains finite. 

By integrating over this region with |q| finite and then taking the limit, we find 

d 3 p 
(2n) 3 2E' 

so that pf(u,p) effectively becomes a delta function. Thus the final result is that the spectral 
density reduces to 



/■q 2 ry r d 3 p 

hmJ Q duj 2 p 2 (uj, q ) = -J -—n F (E)v(3-v) = p sc /8n, (A5) 



lim pf(a;,q) = 5{uj 2 )p sc /8tt . (A6) 



It should be noted here that we did not make any approximation from ( [A4 ) through 
), thus the result is exact. Also note that there arises no ambiguity from the nucleon 
form factor because iq(g = 0) = 1. 
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Figure Captions 



Fig.l: Borel curve for m v (M 2 ) using Tl R . Solid, dashed, dash-dotted lines correspond to 
n B /n = 0, 1.0,2.0 respectively. S (n B ) determined by the Borel stability method at 
each density is also shown in GeV 2 unit. The Borel window is chosen to be 0.41GeV 2 < 
M 2 < 1.30GeV 2 . 

Fig. 2: A schematic illustration of the V — N scattering with s-channel nucleon resonances. 

Fig. 3: (a) Borel curves for m v (M 2 ) at zero density using 11^ with several different values 
of So in GeV 2 unit, (b) Borel curves for m v (M 2 ) at nuclear matter density using U R 
with several different values of S in GeV 2 unit. 
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